A hemoprotein with a zinc-mirror heme site ties heme availability to carbon metabolism in cyanobacteria.

Heme has a critical role in the chemical framework of the cell as an essential protein cofactor and signaling molecule that controls diverse processes and molecular interactions. Using a phylogenomics-based approach and complementary structural techniques, we identify a family of dimeric hemoproteins comprising a domain of unknown function DUF2470. The heme iron is axially coordinated by two zinc-bound histidine residues, forming a distinct two-fold symmetric zinc-histidine-iron-histidine-zinc site. Together with structure-guided in vitro and in vivo experiments, we further demonstrate the existence of a functional link between heme binding by Dri1 (Domain related to iron 1, formerly ssr1698) and post-translational regulation of succinate dehydrogenase in the cyanobacterium Synechocystis, suggesting an iron-dependent regulatory link between photosynthesis and respiration. Given the ubiquity of proteins containing homologous domains and connections to heme metabolism across eukaryotes and prokaryotes, we propose that DRI (Domain Related to Iron; formerly DUF2470) functions at the molecular level as a heme-dependent regulatory domain.

gRNA-SeqRET: a universal tool for targeted and genome-scale gRNA design and sequence extraction for prokaryotes and eukaryotes.

High-throughput genetic screening is frequently employed to rapidly associate gene with phenotype and establish sequence-function relationships. With the advent of CRISPR technology, and the ability to functionally interrogate previously genetically recalcitrant organisms, non-model organisms can be investigated using pooled guide RNA (gRNA) libraries and sequencing-based assays to quantitatively assess fitness of every targeted locus in parallel. To aid the construction of pooled gRNA assemblies, we have developed an in silico design workflow for gRNA selection using the gRNA Sequence Region Extraction Tool (gRNA-SeqRET). Built upon the previously developed CCTop, gRNA-SeqRET enables automated, scalable design of gRNA libraries that target user-specified regions or whole genomes of any prokaryote or eukaryote. Additionally, gRNA-SeqRET automates the bulk extraction of any regions of sequence relative to genes or other features, aiding in the design of homology arms for insertion or deletion constructs. We also assess in silico the application of a designed gRNA library to other closely related genomes and demonstrate that for very closely related organisms Average Nucleotide Identity (ANI) > 95% a large fraction of the library may be of relevance. The gRNA-SeqRET web application pipeline can be accessed at https://grna.jgi.doe.gov. The source code is comprised of freely available software tools and customized Python scripts, and is available at https://bitbucket.org/berkeleylab/grnadesigner/src/master/ under a modified BSD open-source license (https://bitbucket.org/berkeleylab/grnadesigner).

A predicted CRISPR-mediated symbiosis between uncultivated archaea.

CRISPR-Cas systems defend prokaryotic cells from invasive DNA of viruses, plasmids and other mobile genetic elements. Here, we show using metagenomics, metatranscriptomics and single-cell genomics that CRISPR systems of widespread, uncultivated archaea can also target chromosomal DNA of archaeal episymbionts of the DPANN superphylum. Using meta-omics datasets from Crystal Geyser and Horonobe Underground Research Laboratory, we find that CRISPR spacers of the hosts Candidatus Altiarchaeum crystalense and Ca. A. horonobense, respectively, match putative essential genes in their episymbionts' genomes of the genus Ca. Huberiarchaeum and that some of these spacers are expressed in situ. Metabolic interaction modelling also reveals complementation between host-episymbiont systems, on the basis of which we propose that episymbionts are either parasitic or mutualistic depending on the genotype of the host. By expanding our analysis to 7,012 archaeal genomes, we suggest that CRISPR-Cas targeting of genomes associated with symbiotic archaea evolved independently in various archaeal lineages.

Synthase-selected sorting approach identifies a beta-lactone synthase in a nudibranch symbiotic bacterium.

BACKGROUND: Nudibranchs comprise a group of > 6000 marine soft-bodied mollusk species known to use secondary metabolites (natural products) for chemical defense. The full diversity of these metabolites and whether symbiotic microbes are responsible for their synthesis remains unexplored. Another issue in searching for undiscovered natural products is that computational analysis of genomes of uncultured microbes can result in detection of novel biosynthetic gene clusters; however, their in vivo functionality is not guaranteed which limits further exploration of their pharmaceutical or industrial potential. To overcome these challenges, we used a fluorescent pantetheine probe, which produces a fluorescent CoA-analog employed in biosynthesis of secondary metabolites, to label and capture bacterial symbionts actively producing these compounds in the mantle of the nudibranch Doriopsilla fulva. RESULTS: We recovered the genome of Candidatus Doriopsillibacter californiensis from the Ca. Tethybacterales order, an uncultured lineage of sponge symbionts not found in nudibranchs previously. It forms part of the core skin microbiome of D. fulva and is nearly absent in its internal organs. We showed that crude extracts of D. fulva contained secondary metabolites that were consistent with the presence of a beta-lactone encoded in Ca. D. californiensis genome. Beta-lactones represent an underexplored group of secondary metabolites with pharmaceutical potential that have not been reported in nudibranchs previously. CONCLUSIONS: Altogether, this study shows how probe-based, targeted sorting approaches can capture bacterial symbionts producing secondary metabolites in vivo. Video Abstract.

Expression and characterization of spore coat CotH kinases from the cellulosomes of anaerobic fungi (Neocallimastigomycetes).

Anaerobic fungi (Neocallimastigomycetes) found in the guts of herbivores are biomass deconstruction specialists with a remarkable ability to extract sugars from recalcitrant plant material. Anaerobic fungi, as well as many species of anaerobic bacteria, deploy multi-enzyme complexes called cellulosomes, which modularly tether together hydrolytic enzymes, to accelerate biomass hydrolysis. While the majority of genomically encoded cellulosomal genes in Neocallimastigomycetes are biomass degrading enzymes, the second largest family of cellulosomal genes encode spore coat CotH domains, whose contribution to fungal cellulosome and/or cellular function is unknown. Structural bioinformatics of CotH proteins from the anaerobic fungus Piromyces finnis shows anaerobic fungal CotH domains conserve key ATP and Mg2+ binding motifs from bacterial Bacillus CotH proteins known to act as protein kinases. Experimental characterization further demonstrates ATP hydrolysis activity in the presence and absence of substrate from two cellulosomal P. finnis CotH proteins when recombinantly produced in E. coli. These results present foundational evidence for CotH activity in anaerobic fungi and provide a path towards elucidating the functional contribution of this protein family to fungal cellulosome assembly and activity.

Chromosomal integration of complex DNA constructs using CRAGE and CRAGE-Duet systems.

Our recent development of the CRAGE (chassis-independent recombinase-assisted genome engineering) system enables single-step integration of large, complex DNA constructs directly into bacteria genomes across multiple phyla. This protocol describes the details of the experimental design and procedures of CRAGE and extended CRAGE-Duet systems. It also describes a strategy that combines CRISPR with CRAGE, which allows implementation of CRISPR-Cas9, CRISPRa, and CRISPRi in diverse bacteria, overcoming major limitations to broaden the application of CRISPR in non-model bacterial genome engineering. For complete details on the use and execution of this protocol, please refer to Wang et al. (2019), Wang et al. (2020), and Liu et al. (2020).

Protocol for single-cell isolation and genome amplification of environmental microbial eukaryotes for genomic analysis.

We describe environmental microbial eukaryotes (EMEs) sample collection, single-cell isolation, lysis, and genome amplification, followed by the rDNA amplification and OTU screening for recovery of high-quality species-specific genomes for de novo assembly. These protocols are part of our pipeline that also includes bioinformatic methods. The pipeline and its application on a wide range of phyla of different sample complexity are described in our complementary paper. In addition, this protocol describes optimized lysis, genome amplification, and OTU screening steps of the pipeline. For complete details on the use and execution of this protocol, please refer to Ciobanu et al. (2021).

Evolution of trimethoprim/sulfamethoxazole resistance in Shewanella algae from the perspective of comparative genomics and global phylogenic analysis.

OBJECTIVE: Shewanella algae is a zoonotic marine bacterium that causes a variety of infections in immunocompromised patients or those who have been exposed to seawater. The development of trimethoprim/sulfamethoxazole (TMP/SMX) resistance in S. algae are found in human and environment isolates during the past ten years, and thus the treatment options are decreasing. METHODOLOGY: In the study, we conduct a comparative genomic study to identify the resistant mechanism of TMP/SMX-resistance in S. algae. RESULTS: We found the resistance of TMP/SMX in S. algae is associated with the existence of sul1 and dfrA12 within the class 1 integron. The gene cassette dfra12-aadA2-qacEΔ1/sul1 within the class 1 integron is highly conserved. In addition, the class 1 integron and encapsulated sul1 are significantly enriched in Enterobacteriaceae in NCBI and UniProt databases. CONCLUSION: Our study suggests that the horizontal transfer of TMP/SMX resistance via class 1 integron is most frequently occurred within Enterobacteriaceae and has spread to a wide range of sources including soil, poultry, and marine water.

CRAGE-CRISPR facilitates rapid activation of secondary metabolite biosynthetic gene clusters in bacteria.

With the advent of genome sequencing and mining technologies, secondary metabolite biosynthetic gene clusters (BGCs) within bacterial genomes are becoming easier to predict. For subsequent BGC characterization, clustered regularly interspaced short palindromic repeats (CRISPR) has contributed to knocking out target genes and/or modulating their expression; however, CRISPR is limited to strains for which robust genetic tools are available. Here we present a strategy that combines CRISPR with chassis-independent recombinase-assisted genome engineering (CRAGE), which enables CRISPR systems in diverse bacteria. To demonstrate CRAGE-CRISPR, we select 10 polyketide/non-ribosomal peptide BGCs in Photorhabdus luminescens as models and create their deletion and activation mutants. Subsequent loss- and gain-of-function studies confirm 22 secondary metabolites associated with the BGCs, including a metabolite from a previously uncharacterized BGC. These results demonstrate that the CRAGE-CRISPR system is a simple yet powerful approach to rapidly perturb expression of defined BGCs and to profile genotype-phenotype relationships in bacteria.

A single-cell genomics pipeline for environmental microbial eukaryotes.

Single-cell sequencing of environmental microorganisms is an essential component of the microbial ecology toolkit. However, large-scale targeted single-cell sequencing for the whole-genome recovery of uncultivated eukaryotes is lagging. The key challenges are low abundance in environmental communities, large complex genomes, and cell walls that are difficult to break. We describe a pipeline composed of state-of-the art single-cell genomics tools and protocols optimized for poorly studied and uncultivated eukaryotic microorganisms that are found at low abundance. This pipeline consists of seven distinct steps, beginning with sample collection and ending with genome annotation, each equipped with quality review steps to ensure high genome quality at low cost. We tested and evaluated each step on environmental samples and cultures of early-diverging lineages of fungi and Chromista/SAR. We show that genomes produced using this pipeline are almost as good as complete reference genomes for functional and comparative genomics for environmental microbial eukaryotes.

Intracellular pathways for lignin catabolism in white-rot fungi.

Lignin is a biopolymer found in plant cell walls that accounts for 30% of the organic carbon in the biosphere. White-rot fungi (WRF) are considered the most efficient organisms at degrading lignin in nature. While lignin depolymerization by WRF has been extensively studied, the possibility that WRF are able to utilize lignin as a carbon source is still a matter of controversy. Here, we employ 13C-isotope labeling, systems biology approaches, and in vitro enzyme assays to demonstrate that two WRF, Trametes versicolor and Gelatoporia subvermispora, funnel carbon from lignin-derived aromatic compounds into central carbon metabolism via intracellular catabolic pathways. These results provide insights into global carbon cycling in soil ecosystems and furthermore establish a foundation for employing WRF in simultaneous lignin depolymerization and bioconversion to bioproducts-a key step toward enabling a sustainable bioeconomy.

Detection of S83V GyrA mutation in quinolone-resistant Shewanella algae using comparative genomics.

BACKGROUND: Shewanella algae is a zoonotic pathogen that poses a serious health threat to immunocompromised hosts. Treatment of S. algae infections is challenging due to the pathogen's intrinsic resistance to a variety of ß-lactam antibiotics. Therapeutic options have become further limited by the emergence of quinolone-resistant strains. Currently, there are few studies concerning the genetic and molecular mechanisms underlying acquired quinolones resistance in S. algae. qnrA was once proposed as the candidate gene related to quinolones resistance in S. algae. However, recent studies demonstrated qnrA are highly conservative and does not confer resistance to quinolones in S. algae. METHODS: A total of 27 non-duplicated isolates of S. algae strains were examined. MICs of ciprofloxacin were determined using Vitek 2. Whole genome sequencing was performed using MiSeq platform. Comprehensive Antibiotic Resistance Database and ResFinder were used for annotation of quinolones resistance genes. Multiple sequence alignment by EMBOSS Clustal Omega were used to identified mutation in quinolone resistance-determining regions. To investigation of the alteration of protein structure induced by mutation, in silico molecular docking studies was conducted using Accryl Discovery studio visualizer. RESULTS: All S. algae harbored the quinolone-resistance associated genes (qnrA, gyrA, gyrB, parC, and parE) regardless its resistance to ciprofloxacin. Comparison of these genomes identified a nonsynonymous mutation (S83V) in chromosome-encoded gyrase subunits (GyrA) in quinolone-resistant strain. We found this mutation disrupts the water-metal ion bridge, reduces the affinity of the quinolone-enzyme complex for the metal ions and therefore decrease the capability of quinolones to stabilize cleavage complexes. CONCLUSIONS: The study provides insight into the quinolone resistance mechanisms in S. algae, which would be helpful for the evolution of antibiotic resistance in this bacterium.

Bacterial genome editing by coupling Cre-lox and CRISPR-Cas9 systems.

The past decade has been a golden age for microbiology, marked by the discovery of an unprecedented increase in the number of novel bacterial species. Yet gaining biological knowledge of those organisms has not kept pace with sequencing efforts. To unlock this genetic potential there is an urgent need for generic (i.e. non-species specific) genetic toolboxes. Recently, we developed a method, termed chassis-independent recombinase-assisted genome engineering (CRAGE), enabling the integration and expression of large complex gene clusters directly into the chromosomes of diverse bacteria. Here we expand upon this technology by incorporating CRISPR-Cas9 allowing precise genome editing across multiple bacterial species. To do that we have developed a landing pad that carries one wild-type and two mutant lox sites to allow integration of foreign DNA at two locations through Cre-lox recombinase-mediated cassette exchange (RMCE). The first RMCE event is to integrate the Cas9 and the DNA repair protein genes RecET, and the second RMCE event enables the integration of customized sgRNA and a repair template. Following this workflow, we achieved precise genome editing in four different gammaproteobacterial species. We also show that the inserted landing pad and the entire editing machinery can be removed scarlessly after editing. We report here the construction of a single landing pad transposon and demonstrate its functionality across multiple species. The modular design of the landing pad and accessory vectors allows design and assembly of genome editing platforms for other organisms in a similar way. We believe this approach will greatly expand the list of bacteria amenable to genetic manipulation and provides the means to advance our understanding of the microbial world.

Comparative genomics reveals insights into characterization and distribution of quorum sensing-related genes in Shewanella algae from marine environment and clinical sources.

Shewanella algae is not only the most commonly reported species in Shewanella human infections but also capable to inhabit a wide variety of habitats. Although there is evidence that quorum sensing is associated with bacterial adaptation to changing environmental conditions, little is known of the quorum sensing system in S. algae. In this study, we conducted the whole genome sequencing of S. algae strains and applied comparative genomics to reveal the core genome. Genes related to the quorum sensing system were identified by integrated bioinformatics analysis. S. algae harbor genes involved in all three main types of autoinducer systems. This study provides insights into the quorum sensing systems in S. algae, which might be valuable in the future study of cell behavior in S. algae.

CRAGE-Duet Facilitates Modular Assembly of Biological Systems for Studying Plant-Microbe Interactions.

Developing sustainable agricultural practices will require increasing our understanding of plant-microbe interactions. To study these interactions, new genetic tools for manipulating nonmodel microbes will be needed. To help meet this need, we recently reported development of chassis-independent recombinase-assisted genome engineering (CRAGE). CRAGE relies on cassette exchange between two pairs of mutually exclusive lox sites and allows direct, single-step chromosomal integration of large, complex gene constructs into diverse bacterial species. We then extended CRAGE by introducing a third mutually exclusive lox site, creating CRAGE-Duet, which allows modular integration of two constructs. CRAGE-Duet offers advantages over CRAGE, especially when a cumbersome recloning step is required to build single-integration constructs. To demonstrate the utility of CRAGE-Duet, we created a set of strains from the plant-growth-promoting rhizobacterium Pseudomonas simiae WCS417r that expressed various fluorescence marker genes. We visualized these strains simultaneously under a confocal microscope, demonstrating the usefulness of CRAGE-Duet for creating biological systems to study plant-microbe interactions.

An Integrated Computer-Aided Design and Manufacturing Workflow for Synthetic Biology.

Biological computer-aided design and manufacturing (bioCAD/CAM) tools facilitate the design and build processes of engineering biological systems using iterative design-build-test-learn (DBTL) cycles. In this book chapter, we highlight some of the bioCAD/CAM tools developed and used at the US Department of Energy (DOE) Joint Genome Institute (JGI), Joint BioEnergy Institute (JBEI), and Agile BioFoundry (ABF). We demonstrate the use of these bioCAD/CAM tools on a common workflow for designing and building a multigene pathway in a hierarchical fashion. Each tool presented in this book chapter is specifically tailored to support one or more specific steps in a workflow, can be integrated with the others into design and build workflows, and can be deployed at academic, government, or commercial entities.

Bidirectional titration of yeast gene expression using a pooled CRISPR guide RNA approach.

Most classic genetic approaches utilize binary modifications that preclude the identification of key knockdowns for essential genes or other targets that only require moderate modulation. As a complementary approach to these classic genetic methods, we describe a plasmid-based library methodology that affords bidirectional, graded modulation of gene expression enabled by tiling the promoter regions of all 969 genes that comprise the ito977 model of Saccharomyces cerevisiae's metabolic network. When coupled with a CRISPR-dCas9-based modulation and next-generation sequencing, this method affords a library-based, bidirection titration of gene expression across all major metabolic genes. We utilized this approach in two case studies: growth enrichment on alternative sugars, glycerol and galactose, and chemical overproduction of betaxanthins, leading to the identification of unique gene targets. In particular, we identify essential genes and other targets that were missed by classic genetic approaches.

An evolutionarily ancient mechanism for regulation of hemoglobin expression in vertebrate red cells.

The oxygen transport function of hemoglobin (HB) is thought to have arisen ∼500 million years ago, roughly coinciding with the divergence between jawless (Agnatha) and jawed (Gnathostomata) vertebrates. Intriguingly, extant HBs of jawless and jawed vertebrates were shown to have evolved twice, and independently, from different ancestral globin proteins. This raises the question of whether erythroid-specific expression of HB also evolved twice independently. In all jawed vertebrates studied to date, one of the HB gene clusters is linked to the widely expressed NPRL3 gene. Here we show that the nprl3-linked hb locus of a jawless vertebrate, the river lamprey (Lampetra fluviatilis), shares a range of structural and functional properties with the equivalent jawed vertebrate HB locus. Functional analysis demonstrates that an erythroid-specific enhancer is located in intron 7 of lamprey nprl3, which corresponds to the NPRL3 intron 7 MCS-R1 enhancer of jawed vertebrates. Collectively, our findings signify the presence of an nprl3-linked multiglobin gene locus, which contains a remote enhancer that drives globin expression in erythroid cells, before the divergence of jawless and jawed vertebrates. Different globin genes from this ancestral cluster evolved in the current NPRL3-linked HB genes in jawless and jawed vertebrates. This provides an explanation of the enigma of how, in different species, globin genes linked to the same adjacent gene could undergo convergent evolution.

Author Correction: Cryptic inoviruses revealed as pervasive in bacteria and archaea across Earth's biomes.

An amendment to this paper has been published and can be accessed via a link at the top of the paper.

Constructing a yeast to express the largest cellulosome complex on the cell surface.

Cellulosomes, which are multienzyme complexes from anaerobic bacteria, are considered nature's finest cellulolytic machinery. Thus, constructing a cellulosome in an industrial yeast has long been a goal pursued by scientists. However, it remains highly challenging due to the size and complexity of cellulosomal genes. Here, we overcame the difficulties by synthesizing the Clostridium thermocellum scaffoldin gene (CipA) and the anchoring protein gene (OlpB) using advanced synthetic biology techniques. The engineered Kluyveromyces marxianus, a probiotic yeast, secreted a mixture of dockerin-fused fungal cellulases, including an endoglucanase (TrEgIII), exoglucanase (CBHII), ß-glucosidase (NpaBGS), and cellulase boosters (TaLPMO and MtCDH). The confocal microscopy results confirmed the cell-surface display of OlpB-ScGPI and fluorescence-activated cell sorting analysis results revealed that almost 81% of yeast cells displayed OlpB-ScGPI. We have also demonstrated the cellulosome complex formation using purified and crude cellulosomal proteins. Native polyacrylamide gel electrophoresis and mass spectrometric analysis further confirmed the cellulosome complex formation. Our engineered cellulosome can accommodate up to 63 enzymes, whereas the largest engineered cellulosome reported thus far could accommodate only 12 enzymes and was expressed by a plasmid instead of chromosomal integration. Interestingly, CipA 2B9C (with two cellulose binding modules, CBM) released significantly higher quantities of reducing sugars compared with other CipA variants, thus confirming the importance of cohesin numbers and CBM domain on cellulosome complex. The engineered yeast host efficiently degraded cellulosic substrates and released 3.09 g/L and 8.61 g/L of ethanol from avicel and phosphoric acid-swollen cellulose, respectively, which is higher than any previously constructed yeast cellulosome.